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ABSTRACT 

After almost 4 years of operation, the two instruments onboard the Fermi Gamma- 
ray Space Telescope have shown that the number of gamma-ray bursts with high energy 
photon emission above 100 MeV cannot exceed roughly 9% of the total number of 
all such events, at least at the present detection limits. In a recent paper (Zheng et 
al. 2012c), we found that GRBs with photons detected in the Large Area Telescope 
(LAT) have a surprisingly broad distribution with respect to the observed event photon 
number. Extrapolation of our empirical fit to numbers of photons below our previous 
detection limit suggests that the overall rate of such low flux events could be estimated 
by standard image co-adding techniques. In this case, we have taken advantage of the 
excellent angular resolution of the Swift mission to provide accurate reference points for 
79 GRB events which have eluded any previous correlations with high energy photons. 
We find a small but significant signal in the co-added field. Guided by the extrapolated 
power law fit previously obtained for the number distribution of GRBs with higher 
fluxes, the data suggests that only a small fraction of GRBs are sources of high energy 
photons. 

Subject headings: gamma-ray burst: general, gamma rays: stars 

1. Introduction 



A most unexpected feature of gamma-ray bursts (GRBs) became apparent in 1994 with the 
Compton Gamma-Ray Observatory (CGRO) detection of an 18 GeV photon associated with GRB 
940217 (Hurley et al. 1994). About a half-dozen high energy burst events were detected throughout 
the CGRO mission but the actual fraction of the total rate was poorly determined (Catelli et al. 
1998, Dingus 2003). Prior to the launch of the Fermi Gamma-ray Space Telescope in 2008, it was 
possible to anticipate that the LAT instrument would detect more than 200 GRB events per year 
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at energies above 100 MeV (Dingus 2003). Since the GRB spectral energy distribution at lower 
energies has been well characterized by a modified power law with peak fluxes at energies of the 
order of 200 KeV (Kaneko et al. 2006, Goldstein et al. 2012), the existence of photons at energies 
10 times higher is a significant constraint on credible models of the GRB phenomenon (e.g. Band 
et al. 2009). Thus, there had been some anticipation that the Fermi mission might considerably 
enrich our knowledge of this aspect of GRB behavior. Such hopes were dampened by the Band 
paper which was drafted 7 months after the Fermi launch. With three LAT-detected events and an 
extrapolation of BATSE spectral data, Band et al. concluded that the high energy detection rate of 
the LAT would be about one per month. Almost four years of operation has demonstrated that the 
detection rate of energetic GRBs is closer to 9 per year. One interesting discovery that has emerged 
is that short as well as long bursts contribute to the population of > 100 MeV emitting GRBs. To 
date, six short events have been so identified, roughly comparable to their fraction within the total 
GRB population (GRB 080905A, 081024B, 090228A, 090510, 110529A and 111117A). This suggests 
that the high energy photon generation process is a characteristic of relativistic jets, independent 
of the details of their specific progenitors. Given the paucity of information about energetic bursts, 
we have embarked on a program to dig out as much information as possible from the available 
data. This has led to a string of papers describing the discovery of seven high energy GRBs with 
signatures too faint to be detected by more conventional statistical techniques (Akerlof et al. 2010, 
Akerlof et al. 2011, Zheng et al. 2012a, 2012b, 2012c, 2012d). The following work delves further 
into this realm to understand better the extent of the association of high energy events with the 
parent class of all GRBs. 



2. Data selection and analysis 

In our first search of LAT data for faint GRBs (Akerlof et al. 2011), we recognized that 
the precise localizations of bursts (crpsF < 5') detected by Swift, INTEGRAL and AGILE were 
powerful constraints for identifying associated high energy photons in the Fermi LAT. More recently, 
we estimated the event intensity distribution and discovered that within the bounds set by limited 
statistics, the distribution was well fit by the number of LAT photons per unit area raised to a 
constant negative fixed power (Zheng et al. 2012c). To probe this empirical result at intensities 
below those required for reliable single event identification, we realized that the standard field co- 
adding techniques of optical astronomy would work nicely with fields triggered by Swift and similar 
instruments. The GRB events were selected in the interval from June 11, 2008, (Fermi launch 
date) through February 29, 2012 from the Swift catalog!! for which the Fermi LAT boresight angle 
was computed from the Fermi spacecraft attitude datso Since the publicly available LAT data 
stream was restricted to Pass 7 after 2011 August 6, we used Pass 6 prior to the switchover date and 
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Pass 7 thereafter. Events were included only if they passed the maximum LAT acceptance angle 
of 74° and no previous GRB correlations had been reported. Although the LAT acceptance drops 
significantly for boresight angles greater than 65° , it is still large enough for burst detections such as 
reported by the LAT team for GRB 100414A (Takahashi et al. 2010) and GRB 120624B (Vianello 
& Kocevski 2012). Swift events were also ignored if the Fermi spacecraft was inside the South 
Atlantic Anomaly (SAA) or was otherwise inactive. A total of 79 events passed these cuts and are 
listed in Table 1. For reference later in this paper, this sample is called "All-Swift" . Since we were 
interested in investigating possible correlations with the intensity of bursts at lower energies, three 
further selections were made, based on the Fermi GBM fluence. The fluence data was obtained 
from either the GBM Fermitrig catalog^] or GCN messages (Barthelmy 2000). The former took 
precedence whenever available. There were 46 events with such GBM data and this set is named the 
" Swift-GBM" sample. Third and fourth samples, "SWi/iJoright.GBM" and " Swift Aim.GBM" 
consists respectively of the 14 events with fluences greater than 10 fierg/cm 2 within " Swift _GBM" 
and the complementary set of 32 with fluences below this bound. Since " Swift -.GBM" is a subset 
of "AlLSWi/i" and, in turn, "Swi/t_bright_GBM" and " Swift _dim_GBM" are disjoint subsets of 
" Swift -GBM" , these four co-add sets are not statistically independent. Nevertheless, there are 
correlations to the low energy GRB fluxes that makes comparisons useful. It should be noted 
that very few Pass 7 GRB events are included in any of these data sets: All_Swift - 5 out of 79, 
Swift.GBM - 2 out of 46, Swift .bright. GBM - 1 out of 14 and SwiftAim.GBM - 1 out of 32. 
Thus, any systematic biases that might arise from somewhat different primary data processing are 
statistically insignificant. 
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Table 1. List of 79 SW/Uriggered GRBs 
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Table 1 — Continued 
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2.1. Matched filter weight computation and significance 

The heart of our signal detection technique is the use of a matched filter to optimize the 
probability for correctly identifying a GRB in the presence of random backgrounds of gamma-rays. 
The matched filter technique is widely employed for time domain detection of radar and sonar 
signals but can be easily extended to more complex signal structures. In the simplest case, assume 
a uniform background noise, n(i), and a signal waveform, s(t). The detection significance for the 
signal is maximized by integrating a matched filter, f(t), over the effective duration of the signal. 
Variational methods show that f(t) = c-s(t) where c is a constant. For situations with non-constant 
backgrounds, the appropriate choice for f(t) is given by the ratio, c • s{t)/n{t). As described in 
Akerlof et al. 2011, the total weight for each photon was computed by the product of four quantities 
representing the estimated signal to background ratio for energy, position on the sky, time of arrival 
and LAT photon detection class. Finally, a somewhat ad hoc weight sharing factor, £, was applied 
to kill events in which the apparent total event weight was heavily dominated by only one or two 
photons, thus favoring the extreme tails of the background distribution. 

The photon data selection criteria for all events is the same as described in our previous work 
(Akerlof et al. 2010, 2011 and Zheng et al. 2012b, 2012c). Briefly stated, the photon energy is 
restricted to the range from 100 MeV to 300 GeV. A zenith angle cut of less than 105° is applied 
as recommended by the LAT team. A detection time window is defined by the interval from TO 
to T0+47.5s where TO is the nominal trigger time. The photons passing these three criteria are 
then used for further analysis. The matched weight for each photon was calculated by the method 
described in Akerlof et al. 2011, equations 1-5. Note that all classes of photons are included 
together in computing the matched weight score using our previously described methodology. A 
weight-sharing factor, £, was computed via equation 8 and combined to yield the quantity, C^ w i- 
The final result is the sum of the matched weights for all single photons multiplied by the weight- 
sharing factor, It is this value, C, that is propagated to all further analysis. Note that for each 
of the four co-add sets described above, the photons for every matching LAT data set were rotated 
to a common axis and the co-added ensemble was handled as if it were a single GRB event (see 
Figure 1 for a sky map of the Swift-GJiM co-add composite field). 

These various considerations determined how we treated two GRBs with previously reported 
LAT detections, GRB100728A and GRB110625A. GRB 100728A (Abdo et al., 2011) has been 
included in our analysis since although LAT emission was found to be associated with a late time 
X-ray flare, no LAT emission was reported during the GRB prompt phase from to 167 s. In the 
case of GRB 110625A (Tarn et al., 2012), the event was excluded on the basis of the LAT boresight 
angle of 88° when the GRB occurred. 

The first question is whether there is evidence for any significant association of high energy 
photons with these GRB fields. To test this assumption, 1256 sets of random co-added fields were 
constructed with LAT data appropriate for each of the four samples with 79, 46, 14 and 32 GRB 
events respectively. For each of the 1256 random co-added sets, every single constituent field was 
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Table 1 — Continued 
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selected at random time but with center position identical to the corresponding GRB field. These 
random fields also had boresight and zenith angles similar to the corresponding GRB data. The 
distribution of these weights is shown in Figure 2 for the " Swift-GBM" sample of 46. In this case, 
7 random fields out of the 1256 yielded a weight greater than for the actual set associated with 
GRBs. Similar results were obtained for the other samples as well, setting an overall confidence 
level in the neighborhood of 99% that a modest excess of GRB-related photons has been detected. 
A synopsis of these estimates is provided in Table 2. 

The conclusion that a few GRB photons were associated with Swift-triggered events was probed 
in two different ways. It was noted that the fields associated with the 79 GRB events were not 
uniform in intensity as measured by the sum of matched filter weights over the 1256 observations of 
each field taken randomly in time. By comparison of one half of the random data observations with 
the other, a substantial pair-wise correlation was observed for each individual GRB direction. The 
noisiest 8 GRB fields were removed and the previous analysis re-run to see if background astrophys- 
ical sources were a major contamination. The significance levels for the GRB signals compared to 
the random co-added sets remained at the ~99% level for all except the " Sm/i_dim_GBM" sam- 
ple. Secondly, it was noted that the matched filter weight method was moderately sensitive to 
infrequent statistical fluctuations with large weight values. To avoid this possible instability, a 
ranking procedure was invoked that involved two separate steps. First, the 80th percentile weight 
was determined for every set of 1256 random fields associated with each GRB event. Next, a 
score was computed for every random and GRB field based on simply summing the number of two 
highest photon weights above the 80th percentile cut for each field. For the two co-add samples 
with the greatest GBM fluences, " Swift-GBM" and " £w/t_bright_GBM" , there appeared to be a 
significant correlation with the GRB co-add sets surpassing 93% and 97% of the equivalent random 
groups. In summary, we believe that this analysis has found reasonable evidence for a few GRB 
photons within a number of events but this number is remarkably small. 



2.2. Estimating the number of GRB high energy photons 

Since the Swift sample is apparently associated with a non-zero number of GRB photons, the 
next step was to modify the 1256 random fields previously described by adding a fixed number of 
GRB photons taken randomly, with replacement, from the set of 851 photons described in Table 2 
of Zheng et al. 2012c. Since the only parameter of interest is the GRB photon weight, these values 
were chosen by randomly selecting values from the appropriate 851-fold array. This procedure was 
iterated 1000 times for a total of 1256000 co-added fields for each integer number, n, of injected 
GRB photons. The integral distributions of the matched filter weights for the ensemble are shown 
in Figure 2 for pure random background fields (red line) and 1 to 20 injected GRB photons (blue 
lines). To find the most probable photon number for each data set, the integral distribution for 
each curve was estimated at the actual matched filter value. The quoted photon number was 
interpolated from the two values that bracketed 50%. A similar procedure was pursued for the 



Fig. 1. — LAT high energy photon sky map for the photons created by co-adding 46 GRB fields 
triggered by Swift and confirmed by the Fermi GBM. Each field is centered on the GRB direction 
determined by Swift which is indicated by the blue dot. The relative celestial location of each 
photon is shown by a red dot whose size indicates the LAT photon event class (class 1 = small dot, 
class 2 = medium dot, class 3 or 4 = large dot)(Atwood et al. 2009). The photon angular PSFs, 
estimated from the energy, are indicated by the surrounding dotted circles. The green boundary 
circle with a radius of 16.0° provides an angular scale. The plot axes are aligned so that North is 
up and East is to the right. 



±l-cx ranges by finding the crossing points where the integral distributions took on the values of 
0.158655 and 0.841345. The results for each of the four GRB sample sets are listed in Table 2. In 
order of increasingly more stringent GBM fluence requirements, the average number of photons per 
unit area is (0.23, 0.31, 0.53) m -2 , consistent with the high-energy /low-energy correlation shown 
in Figure 5 of Zheng et al. 2012c. (The "5u>i/t_dim_GBM" field has a corresponding flux of 0.30 
photons m -2 .) This fluence trend lends additional support to the evidence of a small but finite 
number of GRB photons in the Swift-triggered co-added fields. 
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Fig. 2. — Cumulative distributions of LAT random fields with to 20 injected GRB photons as a 
function of matched filter weight. The red line shows the distribution of matched weights for 1256 
fields obtained by co-adding 46 LAT images but no additional simulated photons. The blue lines 
show similar distributions for 1256000 simulated fields with 1 to 20 photons added to each. The 
vertical green line indicates the value for the "Swift_GBM" co-add field. 



3. Implications for the number distribution of high energy GRB photons 



In a recently completed paper, we have found that the distribution of GRBs with photons 
> 100 MeV can be characterized by a distribution function that is proportional to the fixed power 
of the number of photons per unit area (Zheng et al. 2012c). This result is shown in Figure 3 
which depicts the distribution of actual GRB events, the parent distribution and the distribution 
modified by the effective area of the LAT detector and analysis event discrimination efficiencies. As 
noted, the inferred number of possible GRB photons in the Swift co-added fields is relatively small. 
Thus for comparison with observation, we can try making the assumption that the GRB number 
distribution can be extrapolated below the detection threshold that originally defined the fit. We 
find that if the event distribution function observed at higher LAT fluxes can be extrapolated to 
lower rates, the inclusion of a fraction of all GRBs is sufficient to accommodate the numbers inferred 
from our analysis. Pursuing such questions will eventually offer a better insight as to whether high 
energy photon emission is a generic feature of all bursts or only a special subset. 



3.1. Computing total photon number from a power law distribution 

We need to find an expression for the total number of photons that can be expected from a 
set of bursts that are distributed in photon number according to a power law. Assume the number 
of GRB events emitting n photons per unit area is given by 

dN v m 
— — = c • n p (1) 

an 

The exponent, p, is approximately -1.8, leading to the conclusion that this mathematical 
form is only valid above some minimum threshold value for n which will be designated as n rn { n . 
Normalizing the integral of dN/dn to unity defines the maximum range over which the power law 
model can operate: 

nrmn = (-^j) 5 * 1 (2) 

For the problem at hand, we need to compute the total number of photons, vsmft, that can be 
expected in the co-added field corresponding to all Swift triggers within the co-added set. Each 
GRB event will be characterized by a specific effective area determined by the boresight angle of the 
GRB with respect to the LAT z-axis. The probability of contributing m photons to the co-added 
field from the i'th GRB event is determined by the Poisson distribution: 

, . f°° (n ■ ai ) m na dN , 

Pi{m)= / i -T dn 3 

./„ . ml dn 



where a, is the LAT effective area appropriate for the i'th GRB event. The total number of 
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photons is then given by 




v swift = ^2 m - pi ( m ") 



(4) 



rn=l i=l 



where 77(777) is the probability that an event with m photons would escape normal detection criteria. 
In practice, this limits the summation over to 8 photons or less. 77(771) was computed by first 
estimating the probability that an event with m photons would have been previously detected by 
the matched filter weight method described in our earlier papers. This corresponds to a matched 
filter weight sum of approximately 15.0. An essentially pure sample of high energy GRB photons 
was obtained from 4 bright GRBs, 080916C, 090510, 090902B and 090926A with 125, 176, 164 
and 177 photons respectively. For each of these burst events, a Monte Carlo program repeatedly 
computed the matched filter weight sum for m randomly selected photons to obtain an average 
probability of exceeding the threshold of 15.0. The independent results for the photon samples from 
the four bright GRBs were substantially similar and the average for each value of m estimates 7/(777), 
the probability of detection by previous searches. 77(777) is computed immediately as 1 — 77' (m). 

Since dN/dn is simply a power of 77, the evaluation of vswift devolves immediately to a sum- 
mation of terms, each of which can be represented by an incomplete gamma function defined by 



Such functions can be found on familiar mathematical packages such as IDL and Mathematica. 
Taken together, equations 4 & 5 link assumptions about the GRB event photon number distribution 
to the actual number of photons that should be observed for a specified ensemble of exposures. In 
this way, the results given in Table 2 will lead to limits for these otherwise inaccessible regions of 
the event photon number distribution. 

3.2. Reconciling the Swift-triggered photons and the extrapolated number 



As noted previously, the lower limit of the range for a power-law distribution with negative 
exponent is bounded by the finite number of events it describes. That condition is satisfied by 
setting the lower bound, n m i n , by equation (2). Thus, the obvious first question is how many 
photons would be expected in the " Swift_GBM" co-add field if all 46 events obeyed this distribution 
over the range from n m in to 00. With n m in defined by equation 2, this can be obtained by application 
of the mathematical relations of equations 4 & 5. For the value of p which best fits the observed 
event distribution shown in Figure 3, the presumption that all GRBs contribute is clearly violated 
at the 97% confidence level as shown in Figure 2 by overpredicting the number of photons, 14, 
instead of the observed value of 6 given in Table 2. 

The obvious next step is to inquire where the lower bound should be set to match the observed 
photon number. To keep the notation unambiguous, this bound will be designated nthresh an d 




(5) 



distribution 
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Fig. 3. — Complement of the cumulative GRB event number distribution, N(n), as a function of 
the LAT high energy photon flux, n. The blue histogram shows the cumulative distribution for the 
LAT-detected events listed in Table 2 of Zheng et al. 2012c. The red curve is the fit to a power-law 
corrected for detection efficiencies. The black lines show the underlying power-law distribution 
in the region of the fit (solid line), extrapolated to satisfy the number of observed sub-threshold 
quanta (dashed line) and a similarly constrained power-law segment with a slightly greater number 
of events with lower fluxes (dotted line). The arrow marked n m j„ indicates where an extrapolation 
of the power law distribution to lower values of the photon flux, n, would reach unity, ie. include 
all GRBs. 
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n thresh > n min- For the most probable values for p and fSwifti the value for nthresh 1S 1-7 pliotons/m 2 
corresponding to 19% of the full GRB sample (n m j n = 0.2 photons/m 2 ). Varying vswift over the 
±1-<T range changes the percentage of the contributing GRBs from 11% to 40%. Keeping vswift at 
its central value but varying p over its ±l-<r range changes these percentages by a factor of roughly 
two. However if p and vswift are simultaneously set to their 1-a upper limits of -1.5 and 10, the 
fraction of GRBs rises to 100%. Admittedly, this extrapolation from present observations relies on 
extrapolation of the photon number distribution to low photon counts. Thus, the analysis should be 
viewed with circumspection but is intended to illustrate the twin constraints of an apparent power- 
law distribution and the limited total number of detected photons associated with a considerable 
number of bursts. 

The abrupt flux cutoff event distribution model described above is an unnatural extreme that 
favors a few events with fluxes near the current lower detection threshold. An alternative is to 
model the differential distribution in the region from zero to the detection threshold as a sum of 
two terms chosen to match the fitted distribution at threshold. With the constraint implied by the 
total of number of events inferred from the co-added analysis, the cumulative distribution plateaus 
at 24% of all GRBs (see dotted curve in Figure 3). The differential distribution for this segment 
is of the form, a + bn q , where a and b are constants, n is the GRB event flux and q is a fractional 
exponent chosen to match the co-added field photon count. This tends to favor a larger number 
of events at lower fluxes. Thus, despite wide differences in model assumptions, the GRB event 
population seems to divide between the ~25% that emit high energy radiation and the 75% that 
don't. We find it exceedingly strange that nature has found a way to limit the number of events with 
high energy flux intensities in the 1 to 4 photons/m 2 range relative to what occurs at higher fluxes. 
Addressing this question might lead to a better understanding of the GRB high energy emission 
process. Finally, the stringent limits for total photon numbers in the various co-added trigger fields 
strongly suggest that there are very few modest intensity bursts that have gone undetected because 
of inadequate signal processing techniques. 

4. Summary 

In this paper, we have found a small but statistically significant number of GRB-associated 
photons for co-added sets of events triggered by the Swift instrument. Although this number has 
large statistical errors, it points to a conclusion that high energy photon emission is relatively 
infrequent and can not be hidden in the somewhat fainter burst events. The necessary initial 
conditions that give rise to this phenomenon are not at all understood. With a larger event sample, 
it may be possible to better classify the population of high energy emitting GRBs and test models 
that explain their behavior. If nothing else, this analysis again demonstrates the importance of 
GRB detection with good angular resolution. 

We thank Timothy McKay for constructive suggestions for this manuscript. This research is 
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Table 2. Summary of Swift-triggered co-added fields 



data set 


number of matched 
co-added filter 
events weight 


estimated photon # 
photon ±l-er 
number range 


effective" 
area 
(m 2 ) 


GBM 6 
fluence 
(//erg cm -2 ) 


AlLSwift 


79 


587.4 


8.52 4.54 - 13.54 


36.23 




Swift.GBM 


46 


370.0 


6.15 3.15 - 10.15 


19.89 


500.7 


Siiri/*_bright_GBM 


14 


111.8 


3.21 1.54 - 5.60 


5.96 


396.7 


Swi /i_dim_GBM 


32 


228.9 


4.33 2.03 - 7.51 


13.93 


104.0 



a Sum of the individual effective areas of all GRBs in the sample. The areas are 
determined from the LAT boresight angles at the burst trigger time. 

b Sum of the GBM fluences of all GRBs in the sample. These num- 
bers are obtained from from the Fermi GBM Burst Catalog, FERMIGBRST 



(http://heasarc.gsfc.nasa.gov/W3Browse/all/fermigbrst.html), and the GCN catalog 



(searchable via GRBlog at http://grblog.org/grblog.php). 



